Quantum entanglement between a nonlinear nanomechanical 
resonator and a microwave field. 



Charles P Meaney, 1 Ross H McKenzie, 2 and G J Milburn 1 

1 Centre for Quantum Computer Technology, 
School of Mathematical and Physical Sciences, 
The University of Queensland, St Lucia, QLD 4072, Australia 
2 Department of Physics, School of Mathematical and Physical Sciences, 
The University of Queensland, St Lucia, QLD 4072, Australia 

Abstract 

We consider a theoretical model for a nonlinear nanomechanical resonator coupled to a super- 
conducting microwave resonator. The nanomechanical resonator is driven parametrically at twice 
its resonance frequency, while the superconducting microwave resonator is driven with two tones 
that differ in frequency by an amount equal to the parametric driving frequency. We show that the 
semi-classical approximation of this system has an interesting fixed point bifurcation structure. In 
the semi-classical dynamics a transition from stable fixed points to limit cycles is observed as one 
moves from positive to negative detuning. We show that signatures of this bifurcation structure are 
also present in the full dissipative quantum system and further show that it leads to mixed state 
entanglement between the nanomechanical resonator and the microwave cavity in the dissipative 
quantum system that is a maximum close to the semi-classical bifurcation. Quantum signatures of 
the semi-classical limit-cycles are presented. 

PACS numbers: 85.85.+j; 42.50.Wk; 82.40.Bj; 84.40.Dc; 85.25.-j 



I. INTRODUCTION. 



The steady states of driven, dissipative quantum systems can manifest entanglement 
between the component physical systems [IHB]. In such cases quantum correlations survive 
despite dissipation, although the resulting states are not necessarily pure. Steady state 
entanglement, in contrast to entanglement through unitary dynamics, does not depend on 
the initial state. For these reasons, the prospect of engineering driven dissipative systems 
so that their steady states exhibit some desirable entanglement is an intriguing prospect [7]. 

In recent years, another aspect of engineered quantum systems has arisen in the 
context of circuit quantum electrodynamics (circuit QED)[S], nano- mechanics [5] and 
optomechanics[TO] . A characteristic feature of the description of these systems is an effective 
quantisation in which quantum theory is used to describe collective macroscopic degrees of 
freedom rather than atomic degrees of freedom. Typical examples include the voltage and 
current in an equivalent circuit description of a a circuit containing superconducting junc- 
tions or a bulk flexural mode of a nanomechanical resonator. This approach works because 
at low temperatures these collective degrees of freedom largely factor out of the microscopic 
degrees of freedom which remain only in so far as a source of dissipation and noise. 

Of particular interest in this paper is an experimental context involving nano-scale res- 
onators coupled to superconducting coplanar waveguides, these devices are made from alu- 
minium on a silicon substrate, and are placed in a dilution refrigerator. There is ongoing 
effort to cool these nanomechanical systems to close to the ground state of a collective bulk 
flexural mode where quantum mechanical phenomena such as entanglement become mani- 
fest. Rapid experimental progress means that such a quantum regime is now accessible [TT]. 
In this paper we propose a particular kind of engineered dissipative quantum system, a non- 
linear nanomechanical resonator coupled to a superconducting microwave resonator, that 
exhibits an entangled dissipative steady state. The correlations implicit in the entangled 
state are enforced by the correlations between fixed points of the corresponding semiclassical 
dynamical system. 

A nanomechanical resonator can form one plate of a capacitor coupling a coherent driving 
field to the cavity field [12, 13J. The nanomechanical element can also exhibit a significant 
Duffing nonlinearity |14j . Such properties allow the investigation of a nonlinear nanome- 
chanical system. The classical model possesses distinct nonlinear phenomena in its steady 
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state, and we seek a signature of these in the quantum system. 

The paper is structured as follows. In section [TT] we introduce the nonlinear nanome- 
chanical system considered in this paper. We give a description in terms of an effective 
Hamiltonian and a Markov master equation describing dissipation of both microwave and 



mechanical modes. In section |III| we present a detailed analysis of the semiclassical dis- 
sipative dynamics of this system in terms of the fixed point stability. In section |IV| we 
consider the quantum version of the model. We numerically determine the steady state of 
the system in the number basis for both resonators, suitably truncated. The steady state 
entanglement between the mechanics and the field can then be calculated in terms of the 
log negativity. From the steady state density matrix in the number basis, we construct the 
marginal Wigner functions for mechanical and cavity degrees of freedom. We show that, 
as the control parameters are varied through the values at which the semiclassical model 
shows bifurcations, the Wigner functions become double peaked with dominant support on 
the semiclassical fixed points. This enables us to elucidate the nature of the entanglement 
in the quantum steady state. Finally in section [V] we summarise our results and suggest 
new directions for further work. 



II. THE DISSIPATIVE CASSINIAN OSCILLATOR MODEL 

We consider a nonlinear nanomechanical resonator coupled to a superconducting mi- 
crowave resonator. The nanomechanical nonlinearity is of the well-established Duffing (quar- 
tic) type [15]. We start by defining the field mode annihilation operator a of the supercon- 
ducting microwave resonator at frequency u c . We also define two field quadrature operators 
x a and y a as 

x a = x (a + a ] ) , a = x a + iy a , 

1 C 1 ) 
Va = -ir { a - flt ) > a ] = x a - iy a . 

Similarly, we define the mechanical mode annihilation operator b at frequency u m . The 

oscillator mass is defined here to be m, and we also define the position and momentum 
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operators q and p as 

as well as two normalised operators Xb and y& as 

Xfe = o ( & + &t ) ' & = x 6 + k/ b , 

1 ( 3 ) 
y b = -\-{b -b ] ) , b ] = x b -iy b . 

The relevant commutation relations for the microwave and mechanical resonators are 

[a,a f ] = /, [x a ,y a ] = i-I, 

\ (4) 
[6, 6 f ] = I, [ac 6 , y 6 ] = i-J, [q, p] = iM , 

where I is the identity. 

The capacitive coupling between the nanomechanical resonator and the cavity field is 
such as to shift the cavity frequency by an amount proportional to the displacement of the 
nanomechanical resonator. As the bare cavity frequency is very high (GHz) compared to 
the observational time scales of interest we can approximate the interaction so that it is 
proportional to the number operator for the cavity field times the mechanical displacement. 
The mechanical resonator thus acts as a phase modulator for the cavity field. Such a 
coupling can also be achieved in an opto-mechanical context where radiation pressure moves 
one mirror of an optical cavity; see [16] for a derivation. 

In our model, the microwave cavity is coherently driven with two tomes, as in [17J; and the 
mechanical oscillator is driven parametrically. Such a parametrically driven oscillator can 
achieve squeezing which is maximal at the threshold for parametric oscillation [TS]. Finally, 
we include in our model of the nanomechanical an anharmonic component. The dominant 
nonlinearity is quartic (in the Hamiltonian, cubic in the force) in the displacement of the 
nanomechanical resonator, and is known as the "Duffing" nonlinearity; a derivation for 
nanomechanical resonators is found in |19j . 

We thus have the Schrodinger picture Hamiltonian describing the coupled system 

H = hw c a^a + hy (e*ae [ ^ + QaV^) + ^q 2 + —p 2 + Xq 4 

i=l 

+ K cos (2u m t) q 2 + HGoa^aq , (5) 
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where: u c is the frequency of the microwave cavity; e\ and €2 are the amplitudes of the 
two linear drives on the microwave cavity; m and u m are the mass and frequency of the 
mechanical oscillator respectively; X is the amplitude of the Duffing nonlinearity of the 
mechanical resonator; K is the amplitude of the parametric pumping of the nonlinear me- 
chanical resonator; and Go is the amplitude of the coupling between the photon number of 
the cavity and the position of the mechanical resonator. A derivation of the form of the 
capacitive coupling between microwave resonator and mechanical oscillator as a^aq is given 
in [13]; a typical achievable value for the coupling strength is G = 27T x 1.16kHznm _1 [20J. 
An equivalent circuit schematic of the system is given in Figure [TJ A physical schematic of 
the microwave-nanomechanical coupling is given in Figure |2j The reader is referred to [21] 
for a scanning electron microscope (SEM) image of such a device. Nanomechanical Duffing 
oscillators have been discussed in [151 E2]; parametric excitation in the nano-mechanical 
context has been discussed in [131 1231 [21] . 

We now move to an interaction picture in a rotating frame with respect to the average of 
the linear drive frequencies in the microwave resonator space, and the mechanical resonator 
frequency in its space. We set UJ2 ~ UJl = uj m and make the rotating wave approximation by 
ignoring terms with frequency 2u m or above. Next, we linearise about the steady state, 
using the same ansatz as Woolley et al [13], and set e\ = ee -1 ^ and e<i = — ee 1 ^ = — where 
eel. Finally, we choose the phase ip — | to obtain the effective Hamiltonian 

U E = huata + (&V 2 ) + K- (b 2 + tf 2 ) - ih 9 - (a + a*) (b - 6 f ) , (6) 

where the notation (• • • )s ^ s use d for a symmetrised product, such that 

(Vot 2 ) = 1 (oV 2 + bb ] bb ] + bb ]2 b + 6 t 6 a 6 t + b ] bb ] b + 6 f V) , (7) 

the cavity detuning uj is 

w = w c , (8) 

and Xi K i an d 9 are the rescaled mechanical Duffing nonlinearity, parametric pumping am- 
plitude, and cavity-mechanical coupling respectively, where 



X=7Z 72 X > K , 9 = \ hi °0- ( 9 ) 

(2mw m ) 2mu m u m V 2mu n 

Classical trajectories from the nanomechanical portion of this Hamiltonian are the ovals 
of Cassini, and the mechanical system is hence sometimes described as the "Cassinian" 




cavity driving 



co planar cavity 

FIG. 1: A schematic of the system under consideration. A nonlinear nanomechanical resonator 
(NR) of frequency oj m is coupled to a superconducting microwave resonator of frequency uj c . The 
nanomechanical oscillator is driven parametrically at frequency 2co m ', and the microwave cavity 
is driven at two tones oj\ and ui2- The annihilation operators for the microwave and nanome- 
chanical modes are a and b, respectively. In terms of these operators, the coupling takes the 
form a) a [b + &'). The superconducting microwave resonator is modelled as a lumped element LC 
circuit, where Lt and Ct are the inductance and capacitance of the tank circuit respectively. 

oscillator; the quantum version has been previously studied by Wielinga et al [25] . A similar 
model has been extensively discussed by Dykmann and his collaborators [26] and also by 
Peano and Thowart |27j. A typical value for the Duffing nonlinearity is x ~ 6.8 x 10~ 4 s _1 

For a realistic description we adopt a dissipative model. We model both the microwave 
cavity resonator and the mechanical resonator as being damped in zero temperature heat 
baths. A zero temperature heat bath for the cavity is certainly justified as the typical mi- 
crowave cavity is at mK temperature and thus very close to zero |28j . The zero temperature 
heat bath for the nanomechanics is not as good an approximation. However, the mean 
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FIG. 2: A physical schematic of the system under consideration. A nonlinear nanomechanical 
oscillator of frequency uj m is capacitively coupled to a superconducting microwave resonator of 
frequency oj c . A theoretical schematic of the system is shown in Figure [T] 

thermal occupation of the bath n^O does not enter the semi-classical equations, and thus 
the semi-classical bifurcation structure studied in section II III will be the correct one. Yet 
n / does affect the quantum steady state and the quantum entanglement calculation. 
For this paper however, we do not study the temperature dependence of the steady state 
entanglement. We thus consider the zero temperature heat bath for the nanomechanics to 
be satisfactory for the purpose of showing the classical-quantum correspondence of sections 



III and IV The amplitude decay for the microwave cavity is p, and for the nanomechanical 
resonator is 7. We then describe the dissipative dynamics with the master equation (with 
weak damping and the rotating wave approximation for the system-environment couplings) 



dp i 



U E ,p + p {Zapa 1 - a ] ap - pa ] a) + 7 (2bpb j - b ] bp - pb ] b) , (10) 
where p is the density matrix of the coupled system. 



III. SEMI-CLASSICAL FIXED POINT STRUCTURE 



A. Reparameterisation 



For the model denned by the master equation (10), we will now find semi-classical steady 
states. We assume that all coupling parameters are non-zero (uj,x,K,g 7^ 0), and that the 
two dissipation parameters are also positive (/1, 7 > 0). 
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First, we define four convenient dimensionless parameters: a renormalised inverse square 
microwave-nanomechanical coupling parameter £; a renormalised parametric drive param- 
eter k'\ a renormalised microwave cavity dissipation parameter //'; and a renormalised 
nanomechanical dissipation parameter 7'. These are defined by 

2 (to 2 + /i 2 ) 



(ii) 



« = - , 

, (J> 
fi = -, 

00 

, 7 
7 = -• 

00 

In terms of our previous assumptions, all of these new parameters are non-zero («', //, 7' 7^ 0), 
and the renormalised coupling is also positive (£ > 0). Additionally, we note that since 
the original dissipation parameters fi and 7 are always positive, the sign of the corre- 
sponding renormalised parameters is equal to the sign of the detuning to; specifically 
sgn (00) = sgn (//) = sgn(7'). In what follows we will write sgn (00) whenever this value 
is required to remind us of its physical origin. We also define one additional parameter 0, 
which we will see determines the locations of the semi-classical fixed point bifurcations, 

= (1 + £«') 2 - (£Y) 2 • (12) 

In the consideration of the semi-classical fixed point structure, it turns out to be most 
useful to consider the re-scaling of the mode operators 

VI (13) 
b ->■ rjb , 



where 

7/ 



and thus the commutation relations become 

to 2 + n 



H ^(l + sgnHv^), (14) 



b,b j 



77 2 



(15) 



This re-parametrisation is a real scaling if 77 > 0, which is true if both 

> and sgn (x) (^sgn (00) + ^J~(p s j > . (16) 
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Performing this re-parametrisation, the master equation ( 10 ) becomes 
dp 2rj 2 uj 



dt 
where 



a'ap — pa/ a) H — — \ 2bpb — b bp 



pPb 



(17) 



% = a) a + 1 + SgnM ^ (b 2 P 2 ) + & ft 2 S t2 
4 V /s 4 



.sgn (u) a/1 + A 4 ' 



(a + a)) (b - fc f ) . (18) 



We now define new quadrature operators for the re-parametrised modes as 

1 



x. 



„ - (a + a ] ) 



j/6 



-i± ft - P 



(19) 



The quantum equations of motion for both quadratures of both oscillators are then obtained 



from the re-parametrised master equation model (17), and are found to be 



d (X a ) 


= ±[ 




dt 






d (Va) 






dt 


" 2-q 2 




d (x b ) 


1 




dt 


" v 





(z ) - p! (y a ) + sgn (w) \Jl+ p' 2 (y b ) 



+ + sgn (w) y/$) « (z 6 2 ft) s ) + <2/6 3 » 



(20) 



d (y b ) _ 1 
dt 2?7 2 

These equations of motion for the expectations are seen to couple to the equations of motion 
of third order moments, which in turn couple to infinite orders. The only case where this 
does not happen would be if 1 + sgn (u) y/4> = 0, and this case is prevented by our conditions 



Zi (Vb)- (l+sgnH yf$) {(x- b 3 ) + ((x b y b 2 ) s )) 



for a valid real re-parametrisation (16). 
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B. Semi-classical model 



To proceed, we define the semi-classical equations by factorising the third order moments, 
and thus obtain a closed system of equations. Specifically, we make the four assumptions 



(x b 3 ) = (x b ) 3 
(x b 2 y b ) s ) = (x b ) 2 (y b ) 
({x b 3 y b 2 ) s ) = (x b ) (y b y 

(y b 3 ) = (y b ) 3 ■ 



(21) 



This is equivalent to saying that the covariances are much less than the corresponding 
product of the means in each case. After factorising the third order moments, four semi- 
classical variables are then defined by 



(.Xa) ' ^ X a , 

(y a ) >->■ y a , 

(x b ) ^ x b , 

(yb) ^ Vb ■ 

These replacements generate the semi-classical equations of motion 



(22) 



dx a 
dt 


2rf 1 


-fjfx a + y a ] , 


dy a 
dt 


" 2rj 2 


-x a - n'y a + sgn 


dx b 


1 


sgn (u) Jl + fi' 2 x 


dt 


~ 2^2 


dyb 
dt 


1 

~ 2^f 


-£tfx h - - 



(23) 



These semi-classical equations of motion (23) have fixed points (%*■ = = ^ = -#■ = 0) 



which are the steady states of the semi-classical dynamical system if they are stable. Thus, 
we must also consider the stability of the fixed points. The defined semi-classical variables 
can be considered as a vector x, such that about a fixed point x we have 



3c x \x a x ao , y a x b x bo , y b y b o\ > 

-^<5x = M<5x, 

dt 



(24) 
(25) 
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where the Jacobian matrix M is 



M = 





-y 


2 












1 


2 


2^ 









| sgn (w) y/l + /J,' 2 




\ sgn (u) y/l + fl' 2 





€7' 
2 


+ (l + sgn (u) y/(f) x b y b 


2 


+ 1+s ^)^ (xl + 3 y 6 2 ) 










£«' 
2 


- 1+sg t^ M + 


57' 
2 


- (l + sgn (w) v^) x 6 i/6 



(26) 

Stability of the fixed point requires all the eigenvalues of the Jacobian to have a real part 
less than or equal to zero |29j. A real part of exactly zero indicates marginal stability in 
that parameter direction, where the fixed point is neither attractive nor repulsive. Real 
parts strictly less than zero are attracting fixed points which draw in nearby regions in 
phase space. In general, stability may depend on more coupling parameter combinations 
than those which define the fixed points. This is indeed the case here, and will be discussed 
explicitly in the next section. 



C. Semi-classical fixed points 

There are three different classes of fixed points of the semi-classical equations of motion 



(23). The fixed points of the semi-classical equations of motion for the microwave resonator 



field, the first two expressions of (23), give the relationship between the nanomechanical 
steady phase space and the microwave resonator field. Explicitly, for all three classes of 
semi-classical fixed points we have, 

x a = — sgn (co>) cos (arctan fjf) y bo , 

(27) 

Va = ~ sgn (w) sin (arctan //) y bo . 
These relations force a correlation between the mechanical and electromagnetic resonator 
at the fixed points. This will be significant when we consider entanglement in the quantum 



steady state in section IV 



Next, to find the scaled momentum of the nanomechanical steady state y bo , and to express 
the nanomechanical components of the steady state in general, it is better to use polar 
coordinates. We define polar coordinates for the nanomechanical phase space r bo and 8 bo in 
the usual way 

Xb = H cos 9 bo , 

(28) 

yb = n sin 9 bo . 



11 



In polar coordinates, the nanomechanical components of the three classes of semi-classical 
fixed points are 

(rj,tan0 6o ) = (0,-), 

1 + -1-£k'-a/0 n 



l + sgn( W )^' H /' ( 29 ) 

1~V$ -l-f/c' + v^ 



1 + sgn (w) a/0' £7' 
The first class of fixed points is the trivial solution at the origin where both quadratures 
of both the microwave and nanomechanical resonators are zero. The fixed point at the 
origin is sometimes, though not always, stable as will be discussed further below. The 
second and third classes of fixed points are both pairs of fixed points at anti-podal positions 
in the nanomechanical phase space. In both cases, the microwave resonator phase space 
components are also a pair of anti-podal fixed points which are separated by a distance 
proportional to the momentum quadrature of the nanomechanical resonator. The difference 
between the second and third classes of fixed points is that one exists as a pair of fixed 
points on the unit circle for all > (this is the second listed pair for positive detuning 
uj > 0, and the third listed pair for negative detuning uj < 0). This class of fixed points is 
typically stable when it exists (below we will be more specific), and we will call these fixed 
points the "unit circle pair" . We can gain some physical insight as to what is going on here 



by recalling the definition of in (12); we see that > corresponds to 

(1 + in'f - (£ 7 ') 2 > . (30) 

This is simple in the case £k! ^> 1 corresponding to a strong parametric drive (for £ fixed). 
Thus, the change from > to < corresponds to the mechanical damping dominating 
the parametric drive. The remaining class of the second and third classes (the third listed 
pair for positive detuning, and the second listed pair for negative detuning) only exists for 
< < 1. At = these fixed points appear (or annihilate for decreasing rather than 
increasing 0) with the unit circle pair; and then head either towards the origin (for positive 
detuning) or out to infinity (for negative detuning) where they then disappear (or appear for 
decreasing rather than increasing 0). This class is always unstable (again more on stability 
below), and we will call these fixed points the "off-circle pair". Physically, < < 1 
corresponds to either weak damping and weak parametric pumping, or, if they are held 
fixed, to £ being small and hence strong coupling. 
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The important observation to make at this point is that the existence and the positions of 
the fixed points depend upon only two parameters, and £7' (since is a function of only 
these two parameters, and sgn (cj) = sgn(£7')). In fact, the existence of the fixed points 
depends only upon the single parameter 0. Additionally, one reason why this particular 
parametrisation was chosen is now quite clear — it places the possibly stable pair of fixed 
points on the unit circle. We plot the radial and angular components of the nanomechanical 
components of the fixed points in figures [3] and [4] respectively. The colours in these figures 
show the stability which is now discussed. 

We also note that a — > —a, b — > — b is a symmetry of the system, and thus pairs of 
anti-podal fixed points are the expected semi-classical result in both factor spaces. 

The stability of the fixed points is not quite as simple. We see that although the existence 
and the positions of fixed points depend upon only two parameters and £7'), the stability 
is dependent on four: and £7'; but also // and £. Thus, a bifurcation or stability 
diagram must be four-dimensional. Nonetheless, we can project upon to the two-dimensional 
space spanned by and £7' to discuss the bifurcation structure somewhat effectively. 
Importantly, to investigate all of the potential switching of stabilities, a large number of 
trials of random sampling in the four-dimensional parameter space was used in addition to 
analytic methods. This numeric method used a logarithmic distribution spanning 7 orders of 
magnitude in each parameter. Where either the stability was analytically shown to be stable 
or unstable, or the numeric method yielded > 99% of trials as either stable or unstable, we 
have indicated the stability of the parameter space as stable or unstable respectively. The 
resulting two-dimensional diagram is plotted in Figure [5] Later, we show an example of a 
mixed region with stable and unstable regions depending on all four parameters in Figure 
6. 

Referring to the "phase diagram" of Figure |5j we discuss the various regions (or "phases" ) 
of the semi-classical dynamical system. The plot of the fixed points themselves on the 
same axes is also helpful, and the reader is encouraged to also refer to Figure [3] during the 
discussion. The first point to note is that the upper half-plane indicates positive detuning 
(00 > 0), and that the lower half-plane indicates negative detuning (u < 0). From the 
diagram it is clear that negative detuning causes widespread instability of the semi-classical 
system. Positive detuning is thus a better candidate for comparison to the full quantum 
steady state of the next section, though both will be compared. 
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FIG. 3: (Color online) Radial components of the nanomechanical oscillator amplitude of the semi- 
classical fixed points. The existence and radial components of the semi-classical fixed points are 
functions of only one parameter <f> = (1 + ^k') 2 — (C^f') 2 , a particular combination of the scaled 
normalised parametric pumping £k' and the scaled normalised nanomechanical dissipation rate £7' 
(compare equation (|29|)). Here, the fixed points are plotted against the two component parameters 
for greater clarity, particularly in regards to stability. There are three different classes of fixed 
points: the origin {r 2 Q = 0) exists for all parameter values; the unit circle pair (r^ = 1) exists 
for (p > 0; and the off-circle pair exists for < <fi < 1. The colours of the plot indicate the 
stability: green indicates stable fixed points; checkered red indicates unstable fixed points; and 
striped blue indicates fixed points that may or may not be stable. The stability is dependent on 
two additional parameters, a normalised microwave cavity dissipation //, and normalised microwave 
cavity-nanomechanical coupling £, as the striped blue region requires. Note that since the sign of 
the microwave cavity detuning sgn (uj) = sgn(£7'), the positive ^'-axis of the diagram represents 
positive detuning and the negative ^'-axis represents negative detuning. It is clear that negative 
detuning leads to semi-classical instability. 

We first introduce some dynamical systems terminology for the bifurcations occurring 
along region boundaries. Briefly, a pitchfork bifurcation consists of one fixed point changing 
stability and giving birth to two new fixed points. The pitchfork bifurcation is supercritical if 
two stable fixed points are born when one stable fixed point loses stability, and is subcritical if 
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FIG. 4: (Color online) Angular components of the nanomechanical oscillator amplitude of the 



semi-classical fixed points (compare equation (29)). The fixed point at origin is not plotted for 
the obvious reason that its angular component is undefined. The angular components of the unit 
circle and off-circle pairs of semi-classical fixed points (the two classes of fixed points other than 
the origin) are functions of two parameters, the scaled normalised parametric pumping £k' and the 
scaled normalised nanomechanical dissipation rate £7'. The angular components of the fixed points 
are plotted against these. The colours of the plot indicate the stability: green indicates stable fixed 
points; and checkered red indicates unstable fixed points. The unit circle pair are the stable fixed 
points in the positive £7' (positive detuning) half-plane. It is seen that for values of the scaled 
normalised parametric pumping £k' > — 1, the unit circle pair tends to be increasingly aligned with 
the y-axis. For values of the scaled normalised parametric pumping £/c' < — 1, the unit circle pair 
tends to be increasingly aligned with the x-axis. Note that the apparent discontinuity in crossing 
the £«/-axis is not real, since to cross means a flipping of the sign £7' and the detuning u (which 
causes a pole in the parameters as the detuning crosses 0), and which also flips the sign of 
To make a clean crossing and keep the sign of the sign of the parametric pumping must be 
flipped together with the detuning. It is best to consider positive and negative detuning (positive 
and negative values of the normalised scaled nanomechanical dissipation rate £7') separately. 



two unstable fixed points are born when one unstable fixed point gains stability. Saddle node 
bifurcations consist of the annihilation of two fixed points, leaving none, or alternatively, 
the creation of two fixed points where there were previously none. Also present are Hopf 
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FIG. 5: (Color online) A two-dimensional projection of the "phase diagram" of the semi-classical 
microwave cavity-nanomechanical system. The qualitatively different regions of semi-classical fixed 
point structures are plotted against the scaled normalised parametric pumping and the scaled 
normalised nanomechanical dissipation rate £7'. Stability depends upon these two parameters and 
also an additional two: the normalised microwave cavity dissipation rate //, and the normalised 
microwave cavity-nanomechanical coupling £. The existence and positions of the fixed points 
depend upon only the first two parameters, making this two-dimensional projection quite useful. 
The colours of the regions indicate the stability: green indicates the region contains some stable 
fixed points; checkered red indicates the region contains only unstable fixed points; and striped blue 
indicates the region may or may not contain a stable fixed point (depending on the additional two 
parameters). Note that since the sign of microwave cavity detuning sgn (w) = sgn (£7'), the top half 
of the diagram represents positive detuning and the bottom half of the diagram represents negative 
detuning. It is clear that negative detuning leads to semi-classical instability. The origin is always a 
semi-classical fixed point. It is stable in regions B and G; and may be stable or unstable, including 
undergoing Hopf bifurcations, in regions A, E, and G. The fixed point pair on the nanomechanical 
unit circle exists in regions C, D, E, F, G, H, I, K, and L; and is stable in regions C, D, E, F, 
and G. The off-circle fixed points exist in regions C, E, G, I, K; but are always unstable. A fuller 
explanation of the semi-classical dynamics is given in the text. The upper dotted line shows the 
location in parameter space of the systems shown in Figures [9 11 and 13 the lower dotted line 



shows the location in parameter space of the systems shown in Figures 10, 12, and 14, The lines 
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are close to the x-axis for ease of quantum numerical simulation as explained in Section IV 



bifurcations of the fixed point at the origin. A Hopf bifurcation consists of a fixed point 
changing stability and giving birth to a limit cycle. A Hopf bifurcation is supercritical if 
a stable fixed point loses stability and gives birth to a stable limit cycle, and is subcritical 
if an unstable fixed point gains stability giving birth to an unstable limit cycle. For more 
information about bifurcations in dynamical systems, see for example [30] . 

For positive detuning, starting at a large negative parametric pumping (and hence a 
negative value of £«;'), as we move from region F through regions E, A, G (or B), and 
C, to D, there is a marked qualitative change in the semi-classical steady state. Ignoring 
for the moment the changes in stability of the origin, there are two fixed points on the 
unit circle which are initially aligned with the x-axis, these gradually rotate before being 
annihilated at the E-A boundary (in a saddle-node bifurcation with an unstable pair) before 
reappearing at the A-G (or B-C) boundary (again in a saddle- node bifurcation with an 
unstable pair) and then gradually rotating towards being aligned with the y-axis. The 
unstable pair with which the annihilation of the unit circle pair occurs, is born and dies 
at the origin in a subcritical pitchfork bifurcation. For negative detuning, starting at a 
large negative parametric pumping (and hence a positive value of since sgn(£ft') = 
sgnwsgn/t), as we move from region L through regions K, J, and I, to H, there is a similar 
change in the semi-classical steady state. However, all of the negative detuned fixed points 
are unstable, so the similarity is only in terms of the fixed point positions. Again ignoring 
the origin, the two fixed points on the unit circle are initially aligned with the x-axis, these 
gradually rotate before being annihilated at the K-J boundary before reappearing at the J-I 
boundary and then gradually rotating towards being aligned with the y-axis. We will seek 
to compare this movement of the fixed points with the changes in the quantum steady state 
phase space in the next section. 

Returning to the changes in stability of the origin, there are potential Hopf bifurcations 
[50] in regions A, E, and G, dependent on the additional two stability parameters, the nor- 
malised microwave cavity dissipation parameter // and the normalised coupling parameter 
£. The Hopf bifurcation causes the stable origin to become unstable, and creates a stable 
limit cycle surrounding it. In order to show this behaviour, we show trajectories of the 
semi-classical system for two different values of // and £ for a particular point in region 
A. Figure [6] shows the stability at a particular point in region A as a function of the two 
additional stability parameters. Figure [7] shows two trajectories on the stable side of the 
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Hopf bifurcation, while Figure [8] shows two trajectories on the unstable side of the Hopf 
bifurcation. We see from figures [7] and [8] that the origin switches from attracting nearby 
trajectories to repelling them, and that an attracting limit cycle is born as the origin loses 
its stability. In the next section we will see what happens to the quantum steady state phase 
space as we move through this semi-classical bifurcation. 




1 2 3 4 5 



FIG. 6: (Color online) A two-dimensional "phase diagram" of the semi-classical microwave 
cavity-nanomechanical system at a particular value of the scaled normalised parametric pump- 
ing £k' = —0.8 and the scaled normalised mechanical dissipation rate £7' = 0.3. This phase 
diagram thus shows the extra two dimensions present at a particular point in region A in the 
two-dimensional projection of Figure [5] of the full four-dimensional phase space. These extra two 
dimensions are the normalised microwave cavity dissipation rate (jf, and the normalised microwave 
cavity-nanomechanics coupling £. The colours of the regions indicate the stability of the fixed point 
at the origin: green indicates the origin is stable; checkered red indicates the origin is unstable. 
For large values of either the mechanical dissipation rate or the normalised coupling, the origin is 
generally stable. 

Lastly, to demonstrate the effect of the unstable fixed points, and in particular to at- 
tempt to get a semi-classical picture for the negative detuning steady state, we investigate 
some trajectories through phase space. Figure [9] shows some trajectories for a linear sweep 
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FIG. 7: (Color online) The nanomechanical components of a trajectory through the semi- classical 
phase space on the stable side of a Hopf bifurcation of the origin, showing a stable attractive 
fixed point at the origin. The parameters for this trajectory are: a scaled normalised parametric 
pumping of = —0.8; a scaled normalised mechanical dissipation rate of £7' = 0.3; a normalised 
microwave cavity dissipation rate of fj,' = 0.3; and a normalised microwave cavity-nanomechanics 
coupling of £ = 2.4. This places the system in region region A of Figure [5j and the green region 
of Figure [6j Qualitatively similar behaviour is found for other system parameter values lying in 
these same two regions; contrast with Figure [8j The trajectory of plot (a) starts close the origin 
at x ao = Ua = %b = Ub = 0.05. The trajectory of plot (b) starts away from the origin at 



1 «o 



Ua = %b = Vb = 0.75. Both trajectories are seen to be attracted to the stable fixed point 



at the origin. Backing out the original couplings necessary for these example parameter values, 
this trajectory is for uj = 3, x = 0.1, n = —1, g = 2.86, \x = 0.9, 7 = 0.38. 



of increasing parametric pumping k (and thus increasing £/«') with all other couplings held 



constant. Figure 10 shows some trajectories for the same linear sweep, with all the same 
couplings, with the exception of a switched sign for the detuning. We note the positive 
detuning trajectories being attracted to the stable semi-classical fixed points, and the neg- 
ative detuning trajectories heading away from the unstable semi-classical fixed points and 
out into surrounding limit cycles. 
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FIG. 8: (Color online) The nanomechanical components of a trajectory through the semi- classical 
phase space on the unstable side of a Hopf bifurcation of the origin, showing a stable attractive 
limit cycle. The parameters for this trajectory are: a scaled normalised parametric pumping of 
£k' = —0.8; a scaled normalised mechanical dissipation rate of £7' = 0.3; a normalised microwave 
cavity dissipation rate of fjf = 0.3; and a normalised microwave cavity-nanomechanics coupling 
of £ = 2.2. This places the system in region region A of Figure [5j and the checkered red region 
of Figure [6j Qualitatively similar behaviour is found for other system parameter values lying in 
these same two regions; contrast with Figure [7j The trajectory of plot (a) starts close the origin 
at x ao = Ua = %b = Ub = 0.05. The trajectory of plot (b) starts away from the origin at 
x a = Ua = Xb = Ub = 0.75. Trajectory (a) shows that the fixed point at the origin has become 
unstable, and trajectories (a) and (b) together that a stable attracting limit cycle has formed. 
Backing out the original couplings necessary for these example parameter values, this trajectory is 
for u = 3, x = 0.1, k = -1.09, g = 2.99, /x = 0.9, 7 = 0.4. 

IV. QUANTUM STEADY STATES. 
A. Quantum steady state phase space. 

In the previous section we showed the fixed point bifurcations of the semi-classical system. 
Here, we investigate whether there is a signature of the semi-classical bifurcations present 
in the full quantum description. We do this by numerically computing the quantum steady 
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FIG. 9: (Color online) Some semi-classical trajectories for positive detuning (u > 0) through a 
linear sweep of increasing parametric pumping. The stable semi-classical fixed points are marked 
with a green circle, and the unstable semi-classical fixed points are marked with a red circle. The 
trajectories start as white/light blue and move towards blue, such that the darker blue is the better 
approximation of the steady state. The axes here are scaled as the original parameters, so that 



they can be compared to the quantum Wigner functions in Figure 11. The trajectories are clearly 
attracted to the stable semi-classical fixed points. Backing out the original couplings necessary 
for the example parameter values plotted for, oj = 2,x = 1,5 = 1.16,// = 0.2,7 = 0.07. The 
parametric pumping k is linearly swept from k = — 8toK = 4in the six plots (along the upper 
dotted line in Figure [5J . There is qualitatively similar behaviour for all systems in the same region 
of the phase space diagram Figure [5] In terms of the parametrisation used in the text, and in the 
phase diagram of Figure [5j we have £7' = 0.2 while £«/ is swept from £k' = —4 to £k' = 2. 



state density operator, and a corresponding Wigner function, in different regions of the semi- 
classical "phase diagram" . It is hoped that by sweeping the coupling parameters through a 
semi-classical bifurcation, there will be a corresponding qualitative change in the quantum 
steady state manifest in the corresponding Wigner function in phase space. This kind of 
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FIG. 10: (Color online) Some semi-classical trajectories for negative detuning (to < 0) through a 
linear sweep of increasing parametric pumping. The stable semi-classical fixed points are marked 
with a green circle, and the unstable semi-classical fixed points are marked with a red circle. The 
trajectories start as white/light blue and move towards blue, such that the darker blue is the better 
approximation of the steady state. The axes here are scaled as the original parameters, so that 



they can be compared to the quantum Wigner functions in Figure 12 . The trajectories are clearly 
repelled by the unstable semi-classical fixed points, out to an outer limit cycle. Backing out the 
original couplings necessary for the example parameter values plotted for, oj = —2,x = 1,5 = 
1.16,// = 0.2,7 = 0.07. The parametric pumping k is linearly swept from k = —8 to k = 4 in the 
six plots (along the lower dotted line in Figure [5]) . There is qualitatively similar behaviour for all 
systems in the same region of the phase space diagram Figure [5j In terms of the parametrisation 
used in the text, and in the phase diagram of Figure [5j we have £7' = —0.2 while £«' is swept from 
£«' = 2 to £«' = -4. 



correspondence principle has proven to be the case for other dissipative nonlinear quantum 
systems [3TT436] . 

To perform the numerical computation of the quantum steady state we use the Quantum 
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Optics MATLAB toolbox [37] . To do this we approximate the infinite basis of each oscillator 
by truncating in the Fock (number) basis at 14 w 3.74 2 . This means that we must choose 
couplings such that the bifurcation takes place sufficiently close to the origin to be accurately 
approximated by the truncation. This is roughly because a coherent state of amplitude a 
has a mean occupation number of |a| 2 . Given the quantum steady state typically (as we 
shall see direct evidence of in this section) has support centred on the semi-classical steady 
state, fixed points far from the origin (high |a|) will produce high occupations and thus 
inaccurate results if we truncate in the Fock (number) basis. 

We hold all couplings equal save the parametric pumping k (and thus £/«'), and move 
along from left to right in the upper half plane of the phase diagram of Figure |5j The 
couplings we choose are the same as those we chose for the semi-classical trajectories in the 
previous section to allow for direct comparison. These were u = ±2, \ — ^-,9 — 1.16,// = 
0.2,7 = 0.07. The parametric pumping k is linearly swept from k — — 8 to k = 4 in the 
six plots that follow for both positive (u > 0) and negative (u < 0) detuning. In terms 
of the parametrisation used in the text, and in the phase diagram of Figure [5j we have 
£7' = —0.2 while £k' is swept from = —4 to £k' = 2. To visualise the microwave 
cavity and nanomechanical factor spaces we choose to look at their phase spaces by using 
the Wigner function of the reduced density matrix for each system. The Wigner function 
is defined as W(x,y) = 4r dz (x — z\p\x + z) e 1_ «~; for more on the Wigner function 
see [TBI |3"8"] . We plot the series of Wigner functions of the nano-mechanical factor space in 



Figure 11 for positive detuning and in Figure 12 for negative detuning. For completeness, 
we also show the corresponding Wigner functions for the microwave cavity factor space in 



figures [13] and [14] for positive and negative detuning respectively. 

For positive detuning, there is a clear signature of the semi-classical bifurcation. For neg- 
ative detuning, the semi-classical fixed points are unstable, though the limit cycle apparent 



in the trajectories of Figure 10 are a signature which appears to be present in the quantum 



steady state phase space as seen in Figure [12} Certainly, the semi-classical approximation 
used in the previous section provides a useful heuristic for understanding the correlations 
inherent in quantum steady states. 
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FIG. 11: (Color online) Density plots of Wigner functions of the nano- mechanical factor space for 
positive detuning (a; > 0) through a linear sweep of increasing parametric pumping. The Wigner 
function W(x, y) is plotted where x and y are the position and momentum quadratures (though 
not scaled as such) of the nanomechanical oscillator. The quantum steady state shows clear signs 
of the semi-classical bifurcations it undergoes. Particular comparison can be made to the semi- 
classical trajectories of Figure |9| Whilst there is evidence of the semi-classical bifurcations, the 
quantum fluctuations wash out the sharp transitions. In particular, the semi-classical bifurcations 
between regions E, A, G, and C of the semi-classical phase diagram [5] are not discenible. However, 
the signature of the overall transition from x-axis aligned density in region F to y-axis aligned 
density in region D is clearly visible. Backing out the original couplings necessary for the example 
parameter values plotted for, oj = 2, % = lj 9 = 1.16, \i = 0.2, 7 = 0.07. The parametric pumping k 
is linearly swept from k = —8 to k = 4 in the six plots (along the upper dotted line in Figure [5]). 
In terms of the parametrisation used in the text, and in the phase diagram of Figure [5j we have 
£7' = 0.2 while £«/ is swept from £«/ = —4 to £k' = 2. 



B. Quantum steady state entanglement. 

We wish to determine the steady state quantum entanglement between the superconduct- 
ing microwave resonator and the nonlinear nanomechanical resonator. The complete system 
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FIG. 12: (Color online) Density plots of Wigner functions of the nano-mechanical factor space for 
negative detuning (cu < 0) through a linear sweep of increasing parametric pumping. The Wigner 
function W(x,y) is plotted where x and y are the position and momentum quadratures (though not 
scaled as such) of the nanomechanical oscillator. Particular comparison can be made to the semi- 



classical trajectories of Figure 10 The major and minor axes of symmetry of the phase-diffused 



rings are seen to be the same as the semi-classical limit cycles. Backing out the original couplings 
necessary for the example parameter values plotted for, u) = —2, x = 1,3 = 1.16, [i = 0.2, 7 = 0.07. 
The parametric pumping k is linearly swept from n = — 8 to k = 4 in the six plots (along the 
lower dotted line in Figure [5]). In terms of the parametrisation used in the text, and in the phase 
diagram of Figure [5j we have £7' = —0.2 while £«' is swept from = 2 to £k' = —4. 



consists of three parts: the superconducting microwave resonator; the nanomechanical res- 
onator; and the environment. A natural measure of such entanglement is the log negativity 
[39] . With access to the truncated steady state density matrix, we can assume the truncated 
matrix is a sufficiently good approximation and directly compute the log negativity as 
the base 2 logarithm of the trace norm of the partial transpose of the bipartite density 
matrix, 

E N (p) = log 2 \\p TA \\ 1 . (31) 
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FIG. 13: (Color online) Density plots of Wigner functions of the superconducting microwave cavity 
factor space for positive detuning (u > 0) through a linear sweep of increasing parametric pumping. 
The Wigner function W{x, y) is plotted where x and y are two quadratures operators of the 
superconducting microwave cavity field. The quantum steady state shows clear signs of the semi- 
classical bifurcations it undergoes. Backing out the original couplings necessary for the example 
parameter values plotted for, oj = 2, x = 1, 9 = 1.16, fi = 0.2, 7 = 0.07. The parametric pumping k 
is linearly swept from k = —8 to k = 4 in the six plots (along the upper dotted line in Figure [5]). 
In terms of the parametrisation used in the text, and in the phase diagram of Figure [5j we have 
£7' = 0.2 while £«' is swept from £k' = —4 to £«/ = 2. 



We compute the log negativity for each quantum steady state in a linear sweep of the 
parametric pumping k whilst holding other parameters constant at the same values used 
to investigate the quantum steady state phase space in the previous section. We note that 
the entanglement entropy of similar nonlinear systems has been examined, though with the 
quartic nonlinearity treated perturbation [4"D~] . 

The quantum entanglement for the positive and negative detuning cases are plotted in 



figures 15 and 16 respectively. In what we consider the more meaningful "phase transition", 



the bifurcation traversals in the positive detuning case of Figure 15, we see that the en- 
tanglement is peaked slightly to one side of the bifurcation centre. Generally speaking, a 
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FIG. 14: (Color online) Density plots of Wigner functions of the superconducting microwave cav- 
ity factor space for negative detuning (uj < 0) through a linear sweep of increasing parametric 
pumping. The Wigner function W(x, y) is plotted where x and y are two quadratures operators of 
the superconducting microwave cavity field. Backing out the original couplings necessary for the 
example parameter values plotted for, u = —2,x = 1)9 = 116, /U = 0.2,7 = 0.07. The parametric 
pumping k is linearly swept from k = — 8 to k = 4 in the six plots (along the lower dotted line in 
Figure [5]). In terms of the parametrisation used in the text, and in the phase diagram of Figure [5j 
we have £7' = —0.2 while £«' is swept from = 2 to £«/ = —4. 

system has more entanglement at criticality; however the exact peak can depend upon the 
entanglement measure. For example, using concurrence to measure entanglement for the 
transverse Ising model shows a peak just off criticality [JT]; while Vidal et al have shown 
that entanglement entropy scaling with block size shows that for a sufficiently large, fixed 
block size, the entanglement entropy peak is always at criticality [4"2] . The fact that the 
position of the quantum entanglement peak is very close to but not coincident with the 
parameter value at which a phase transition occurs is not uncommon in ID physical systems 
[UI16]. This is typically a finite size effect (i.e for systems that are not in the thermodynamic 
limit). However, for other systems, the entanglement peak is coincident with the parameter 
values at which a phase transition occurs: this has been shown in the Dicke model |43] and 
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Lipkin model [44] in particular. See also the model in [45] . 




FIG. 15: (Color online) Quantum entanglement, as measured by the logarithmic negativity, for 
positive detuning (u > 0) through a linear sweep of increasing parametric pumping. The four 
vertical lines show the positions of the semi-classical bifurcations: the first vertical line (the dotted 
red line) is where the two unstable off-circle fixed points appear at the origin (the F-E boundary); 
the second vertical line (the solid red line) is where just created unstable fixed points annihilate the 
unit circle pair (the E- A boundary) ; the third vertical line (the solid blue line) is where the two pairs 
of fixed points re-appear (the A-G boundary); and the fourth vertical line (the dotted blue line) is 
where the off-circle pair vanishes at the origin (the C-D boundary). The quantum entanglement is 
peaked just to one side of the semi-classical bifurcation centre. Backing out the original couplings 
necessary for the example parameter values plotted for, to = 2, x = 1)5 = 1.16, /i = 0.2,7 = 0.07. 
The parametric pumping k is linearly swept from k = —8 to K = 4. In terms of the parametrisation 
used in the text, and in the phase diagram of Figure [5j we have £7' = 0.2 while £k' is swept from 
= -4 to f/s' = 2. 



28 




FIG. 16: (Color online) Quantum entanglement, as measured by the logarithmic negativity, for 
negative detuning (u < 0) through a linear sweep of increasing parametric pumping. The four 
vertical lines show the positions of the semi-classical bifurcations (from right to left): the first 
vertical line (the dotted red line) is where the two off-circle fixed points appear (the L-K boundary) ; 
the second vertical line (the solid red line) is where just created fixed points annihilate the unit 
circle pair (the K-J boundary); the third vertical line (the solid blue line) is where the two pairs of 
fixed points re-appear (the J-I boundary); and the fourth vertical line (the dotted blue line) is where 
the off-circle pair vanishes (the I-H boundary). The quantum entanglement for negative detuning 
is remarkably flat. Backing out the original couplings necessary for the example parameter values 
plotted for, ui = —2,x = 1,5 = 1.16,// = 0.2,7 = 0.07. The parametric pumping k is linearly 
swept from k = —8 to k = 4. In terms of the parametrisation used in the text, and in the phase 
diagram of Figure [5j we have £7' = —0.2 while is swept from ^k' = 2 to = —4. 
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We also calculate the purity of the quantum steady state (the trace of the square of 
the reduced density matrix at its steady state, Tr{p 2 }). We compute the purity for each 
quantum steady state in a linear sweep of the parametric pumping k whilst holding other 
parameters constant at the same values used to investigate the quantum steady state phase 
space in the previous section. 



The purity for the positive and negative detuning cases are plotted in figures 17 and 18 



respectively. We see that the purity is reduced when the separation of the Wigner function 
density from the origin increases. We interpret this as being consistent with the well-known 
notion that as physical separation of two quantum objects increases, the decoherence 
is greater, and they will become more classical (less pure). 

As mentioned, for ease of numerical simulation, we have picked parameters such that the 
steady state phase space density over our two sweeps of parameter values remains largely 
close to the origin in both the microwave and nano-mechanical spaces. That this means 
our sharp semi-classical bifurcations are more significantly blurred by quantum fluctuations 
is one reason we offer for the relatively small variation of both the entanglement and the 
purity over our chosen sweep. 



V. CONCLUSION. 



In this paper we detailed the semi-classical and quantum steady state structure of a 
particular dissipative nanomechanical system. In particular we observed that the semi- 
classical model contains a rich bifurcation structure, and that the remains of this structure 
are still visible in the full quantum mechanical steady state. The steady state quantum 
entanglement was found to be a maximum just to one side of the centre of the semi-classical 
bifurcation. 

Specifically, the semi-classical fixed point structure contains not just one bifurcation, but a 
series of bifurcations, as a parameter is varied. The semi-classical fixed points' existence and 
locations were shown to be dependent on just two dimensionless parameter combinations, 
while their stability was shown to be dependent on these two together with an additional 
two dimensionless parameters. Semi-classically there were also periodic steady states for 
negative values of the detuning parameter. 

The numerically calculated quantum steady states were shown to have clear signatures of 
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FIG. 17: (Color online) Purity of the steady state density matrix, for positive detuning (oj > 0) 
through a linear sweep of increasing parametric pumping. The four vertical lines show the positions 
of the semi-classical bifurcations: the first vertical line (the dotted red line) is where the two 
unstable off-circle fixed points appear at the origin (the F-E boundary); the second vertical line 
(the solid red line) is where just created unstable fixed points annihilate the unit circle pair (the 
E-A boundary); the third vertical line (the solid blue line) is where the two pairs of fixed points 
re-appear (the A-G boundary); and the fourth vertical line (the dotted blue line) is where the 
off-circle pair vanishes at the origin (the C-D boundary). The purity is reduced as the separation 
of the Wigner function density increases. Backing out the original couplings necessary for the 
example parameter values plotted for, ui = 2,x = 1,5 = 1.16, fj, = 0.2,7 = 0.07. The parametric 
pumping k is linearly swept from k = —8 to k = 4. In terms of the parametrisation used in the 
text, and in the phase diagram of Figure [5j we have Q~f = 0.2 while £/«' is swept from £k' = —4 to 
£«' = 2. 
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FIG. 18: (Color online) Purity of the steady state density matrix, for negative detuning (u < 0) 
through a linear sweep of increasing parametric pumping. The four vertical lines show the positions 
of the semi-classical bifurcations (from right to left): the first vertical line (the dotted red line) 
is where the two off-circle fixed points appear (the L-K boundary); the second vertical line (the 
solid red line) is where just created fixed points annihilate the unit circle pair (the K-J boundary); 
the third vertical line (the solid blue line) is where the two pairs of fixed points re-appear (the J-I 
boundary); and the fourth vertical line (the dotted blue line) is where the off-circle pair vanishes 
(the I-H boundary). The purity for negative detuning is remarkably flat. Backing out the original 
couplings necessary for the example parameter values plotted for, uj = —2,x = 1><7 = 1.16,// = 
0.2, 7 = 0.07. The parametric pumping k is linearly swept from k = —8 to k = 4. In terms of the 
parametrisation used in the text, and in the phase diagram of Figure [5j we have £7' = —0.2 while 
£k' is swept from £k' = 2 to = —4. 
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these semi-classical steady state bifurcations. Specifically, the Wigner function representa- 
tion of the quantum phase space was seen to have support on the semi-classical fixed points. 
In addition, where the semi-classical model had no stable fixed points, but instead had a 
periodic steady state, the Wigner function was seen to have support all around the limit 
cycle. This is consistent with the quantum phase space being completely phase-diffused 
around the limit cycle in the quantum steady state. 

The "Cassinian" oscillator is thus an example of a correspondence principle between 
classical dynamics and quantum steady states. This principle is that investigation of the 
dynamics of the relevant semi-classical model gives significant predictive power for the steady 
state behaviour of the full quantum dissipative system. 

Experimentally, variation of the parametric pumping k, whilst holding all other param- 
eters constant and for a positively-detuning cavity drive u > 0, allows tuning through the 
semi-classical saddle-node bifurcations (from regions F to D in Figure [5]). Semi-classically, 
this means the fixed points switch from being aligned along the x-quadrature (for negative 
k) to being aligned the y-quadrature (for positive k). The calculated quantum steady-state 



phase space of section IV indicates that read-out of the microwave cavity steady state should 
see a splitting in the x quadrature as the semi-classical bifurcation is crossed in accordance 



with Figure 13 Such a splitting becomes sharper the closer in parameter space one is 
to the x-axis in the "phase diagram" of Figure [5] (where a smaller change in parametric 
pumping is required to transit regions E, A, G, B, and C). Experimentally, this means that 
the bifurcation is always achievable; however, the stronger the coupling g, the sharper the 
transition. 

Currently achievable experimental values for such a nano-mechanical system are: a mi- 
crowave cavity of frequency u c = 2n x 5GHz and a linewidth /i = 2tt x 490kHz; a nano- 
mechanical resonator of mass m = 2pg, a frequency under tension u m = 2tc x 2.3MHz, and 
a linewidth 7 = 2tt x 20Hz; and a microwave-mechanical coupling Go = 2i\ x 1.16kHznm _1 
[20] , which yields a linearised g that can be increased by stronger driving up to a maximum 
circulating current of 1/iW p2]. In fact the higher couplings Go = 2n x 6.4kHznm~ 1 |12j 
and even Go = 2n x 32kHznm _1 [50] have been achieved. We believe these values to allow 
the quantum signatures of the just mentioned semi-classical saddle-node bifurcations to be 
seen. 

For negative detuning, and for some parameter regimes (some points in the meta-stable 
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region A of Figure [5]), the nano- mechanics classically decays to periodic steady states, and 



the quantum system become completely phase-diffused, as in Figure 12 However, if read-out 



is to be via the microwave cavity, observation of these bifurcations will remain a challenge, 



since the phase-diffused rings are flattened as in Figure 14 This is essentially because of 



some remnant of the correlation observed in (27) 



In conclusion, even if the nano-mechanical experiment does not achieve a very large 
coupling, there are still semi-classical bifurcations which will be transited by varying the 
parametric pumping power. The overall transition from region F to region D of Figure [5] 
will leave a visible quantum signature which is the separation of steady state phase space 
density of the superconducting microwave cavity along one quadrature, as shown in Figure 
[131 The detailed series of semi-classical bifurcations crossed in this transition will however be 
washed out. A stronger coupling will sharpen the quantum signatures of these semi-classical 
bifurcations. 
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